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Ultracold atoms in an ultrahigh-finesse optical cavity are a powerful platform to produce spin 
squeezing since photon of cavity mode can induce nonlinear spin-spin interaction and thus generate a 
one-axis twisting Hamiltonian Ho at = qJx , whose corresponding maximal squeezing factor scales as 
N ~ 2 / 3 , where N is the atomic number. On the contrary, for the other two-axis twisting Hamiltonian 
Htat = q(Jx — Jy), the maximal squeezing factor scales as iV -1 , approaching the Heisenberg limit. 

In this paper, inspired by recent experiments of cavity-assisted Raman transitions, we propose a 
scheme, in which an ensemble of ultracold six-level atoms interacts with two quantized cavity fields 
and two pairs of Raman lasers, to realize a tunable two-axis spin Hamiltonian H = q(Jx + xJy) + 
LJoJ z . For proper parameters, the above one- and two- axis twisting Hamiltonians are recovered, and 
the scaling of N~ x of the maximal squeezing factor can occur naturally. On the other hand, in the 
two-axis twisting Hamiltonian, spin squeezing is usually reduced when increasing the effective atomic 
resonant frequency cco- Surprisingly, we find that by combined with the dimensionless parameter 
x(> — 1), the effective atomic resonant frequency ujq can enhance spin squeezing largely. These 
results are benefit for achieving the required spin squeezing in experiments. 

PACS numbers: 42.50.Dv, 42.50.Pq 


I. INTRODUCTION 


Squeezed spin states, which were firstly introduced by 
Kitagawa and Ueda |l[, are quantum correlated states 
with reduced fluctuations in one of the collective spin 
components [2, 3]. Such states have attracted consider¬ 
able interest because they not only play significant roles 
in investigating many-body entanglement [3,0, but also 
have important applications in atom interferometers and 
high-precision atom clocks 00- In general, there are 
two methods to produce spin squeezing. One is based on 
a one-axis twisting Hamiltonian 77 oat = qJ x i where q 
is the nonlinear spin-spin interaction strength and J x is 
the collective spin operator in the x direction. When 
the initial state is prepared as \J Z = —j) for q > 0 
(I J z ~j) for g < 0), the maximal squeezing factor for 
this one-axis twisting Hamiltonian scales as TV -2 / 3 0, 
where N is the total atomic number and J = N/2. On 
the contrary, for the other two-axis twisting Hamiltonian 
Htat = q{J x — Jy), the maximal squeezing factor scales 
as TV -1 with the same initial state [1]. Since the scal¬ 
ing of TV ” * 1 approaches the Heisenberg limit, implement¬ 
ing the two-axis twisting Hamiltonian in current exper¬ 
imental setups is very important and necessary [9l-lllj. 
A proposing scheme is to transform the one-axis twist¬ 
ing Hamiltonian into the two-axis twisting Hamiltonian 
by applying pulse sequences or continuous driving in the 
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two-component Bose-Einstein condensates min. How¬ 
ever, the experimental realization of such two-axis twist¬ 
ing Hamiltonian is still challenging. 

Ultracold atoms in an ultrahigh-finesse optical cavity 
are also a powerful platform to produce spin squeezing 
since photon of cavity mode can induce nonlinear spin- 
spin interaction and thus generate the required one-axis 
twisting Hamiltonian j 1 C>l l23j . Recently, multi-mode cav¬ 
ities [24[ have attracted much attention both experimen¬ 
tally and theoretically [2bl l2fij . On one hand, these setups 
can be used to explore novel physics, such as the spin- 
orbit-induced anomalous Hall effect IS El , the crys¬ 
tallization and frustration (32|, |33|, the spin glass [ffil 
[37], and the gapless Nambu-Goldstone-type mode with¬ 
out rotating-wave approximation [38|. Moreover, two¬ 
mode field squeezing |39| and unconditional preparation 
of a two-mode squeezed state of effective bosonic modes 
0 have also been achieved by introducing two cavities. 
In the present paper, inspired by recent experiments of 
cavity-assisted Raman transitions (4l|, |42| , we mainly re¬ 
alize a generalized two-axis spin Hamiltonian by two cav¬ 
ities to enhance spin squeezing largely. 

When an ensemble of ultracold six-level atoms inter¬ 
acts with two quantized cavity fields and two pairs of 
Raman lasers, we first realize a two-mode Dicke model. 
In the dispersive regime, we obtain a generalized two-axis 
spin Hamiltonian H = q(J% + %J 2 ) + ujqJ z , where % is 
a dimensionless parameter and uoq is an effective atomic 
resonant frequency. This realized Hamiltonian has a dis¬ 
tinct property that the interaction strength g, the dimen¬ 
sionless parameter x> and the effective atomic resonant 
frequency ujq can be tuned independently. For reason- 
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able parameters, the one- and two-axis twisting Hamil¬ 
tonians are recovered. Numerical results reveal that for 
the standard two-axis twisting Hamiltonian (% = — 1 and 
cjo = 0) 5 the corresponding maximal squeezing factor 
scales as TV -1 , as expected. On the other hand, in the 
two-axis twisting Hamiltonian Ftat, spin squeezing is 
usually reduced when increasing the effective atomic res¬ 
onant frequency uj o- Surprisingly, we find that by com¬ 
bined with the dimensionless parameter y(> —1), the 
effective atomic resonant frequency c^o can enhance spin 
squeezing largely. These results are benefit for achieving 
the required spin squeezing in experiments. 

This paper is organized as follows. Section II is de¬ 
voted to realizing the generalized two-axis spin Hamilto¬ 
nian with independently-tunable parameters. Section III 
is devoted to introducing the spin squeezing factor. Sec¬ 
tion IV is devoted to numerically investigating the max¬ 
imal squeezing factor since for the two-axis spin Hamil¬ 
tonian, analytical results are very hard to be obtained 
[H. The parts of Discussions and Conclusions are given 
in sections V and VI. 


II. MODEL AND HAMILTONIAN 
A. Proposed experimental setup 

Motivated by recent experiments of cavity-assisted Ra¬ 
man transitions (4l|, |42j, here we propose a scheme, in 
which an ensemble of ultracold six-level atoms interacts 
with two quantized cavity fields and two pairs of Raman 
lasers [see Fig. [TJa)] , to realize a generalized two-axis spin 
Hamiltonian with independently-tunable parameters. As 
shown in Fig.[]Jb), six levels consist of two stable ground 
states (|0> and |1)) and four excited states (|ri), |ri), |si), 
and 1 ^ 2 )). Two independent photon modes, whose cre¬ 
ation and annihilation operators are a) (b^) and a (6), 
mediate the |0) <—> \si) and |1) «—» \ri) (|0) <—» |s2) 
and |1) <—y |r* 2 )) transitions (red and blue solid lines) 
with atom-photon coupling strengths g Sl (g S2 ) an d g ri 
(g r2 ), respectively. Whereas two pairs of Raman lasers 
govern the other transitions {|0) <—y |ri), |1) ^y \ s i)} 
and {|0) <—y |r* 2 ), |1) ^ \ s 2 }} (red and blue dashed 
lines) with Rabi frequencies {D ri , D Sl } and {D r 2 , H S2 }, 
respectively. A ri 2 and A Sl 2 are the detunings from the 
excited states. 



(a) 



(b) 


FIG. 1: (Color online), (a) Proposed experiment ally-feasible 
setup that an ensemble of ultracold six-level atoms interacts 
with two quantized cavity fields and two pairs of Raman 
lasers, (b) Energy-level structure and their transitions in¬ 
duced by two photon modes and two pairs of Raman lasers. 
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B. Total Hamiltonian 

The total Hamiltonian illustrated in Fig.[l]can be writ¬ 
ten as 

H t = He + Ha + H h (1) 

where 

He = uJaa'a + cob^b, ( 2 ) 


N 

Hi = [&»i l s i>j (°lj a + 9 ri \ r i)j (l|j a (4) 

3 = 1 

+ 9 s 2 \* 2 )j (Olj b + g r2 \r 2 )j ( 1 ^ b + H.c. . 

In the Hamiltonians ©-®, uo ri 2 , cj S 1j2 , and uj\ are the 
atomic frequencies, uj\ 312 (ip s12 ) an d ^ln 2 {Vr 12 ) are 
the frequencies (phases) of Raman lasers, respectively, 
and H.c. denotes the Hermitian conjugate. 
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C. Two-mode Dicke model 


By means of the Hamiltonian we first realize a 
two-mode Dicke model with independently-tunable pa¬ 
rameters. In the interaction picture with respect to the 
free Hamiltonian 

N 

H 0 = u^ffl t a + ^ 6 t & + y^ (ay ki), <^i1 7 (5) 

j=i v 

fa)j (ri\j + w 7 2 | s 2 )j {s2\j 
+^v 2 fa)j fa\j +wl |l)j- <i|j) , 
where s = (w 0ri + 0 J\ Sl )/ 2 , uT b = w)f 2 = (wor 2 + 

^ls2 ) /^ 5 ^1 = (^Ori <^lsi)/2 — (^0r2 kq^)/^ ^r\ = 

cJi + a/^, and = cJi + cJJ, the Hamiltonian 0 can be 
rewritten as 


H = 


N 

A a aXa + A fc & t 6 + ^ |A n Iri),. (n^ 


(6) 


3=1 


+ A r 2 |^2 ) j <^*2 \j + A S1 |si)j (si|,- + A S2 |s 2 )j (S 2 1j 
+ a i |l)j (l|j + 2 fi-i l s i)j (l|j e * Vsi 
+^s 2 l s 2),- (1|,- e* Vs2 + On |f*i) ■ 


p™Pr 


+O r2 fa)j (0|j e* Vr2 + H.c. + g Sl |si) 
+9n Vi)j <l|j a + g S2 \s 2 )j (Olj b 
+9r 2 fa)j (l|j b + H.c.] } , 


where A a = cda - a > a , A b = co b A Sl = cj Si - u 3l , 
A S2 - uo S2 — u S2 , and Ai = uj% — uj\. 

In the large-detuning limit, i.e., |A ri2jSl2 | {fi r i >2 , 
D Sl)2 , ^ 1,25 .q .si oi , all excited states can be eliminated 
adiabatically [43|, and an effective Hamiltonian is ob¬ 
tained by 


H = cuAC^a ± oJB^b ± (JqJ z ± ga^aJ z (7) 

+ [(A ri ae"^i + A r2 6e“^) J_ 

+ (A ai ae-^-i + \ S2 be-^^) J+ + H.c.] , 

where j? = 3 ^/A n + 5r 2 / A r 2 - 9 2 s J^ Sl ~ 9s 2 / A s 2 , 
•4 _= |l)j <0|j, J- = |0},<1|,-, and J z = 

J2jZi(\l)j (°lj ~ l°)j (01,0/2 are the collective spin oper- 
ators, 


u, -a + 1 f n -- +“•• (S) 

"°- Al + 3 aT + aT _ aT“ aT 1 ’ (8) 
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FIG. 2: (Color online). Numerical plots of the time- 
dependent spin squeezing factors £§(£) of both the Hamiltoni¬ 
ans m and (fT6|). with the same initial states | J z = —j). In 
(a) and (b), uja/cjo — 200, ujb/cjo — d=200 [“+” for (a) and 
” for (b)], Ai/c^o — 2, A 2 /CJ 0 = 1, and N = 15. In (c) and 
(d), <jJa/ojo — 200, ujb/ojo — ±200 [“±” for (c) and ” for 
(d)], Ai/c^o = 1, A 2 /CJ 0 1 2, and N — 20. 


are the effective atomic resonant frequency and the ef¬ 
fective frequencies of two photon modes a and 6 , respec¬ 
tively, and 


At> ■ - 


]L Qn^rj \ 

2 A ri 5 


\ 9sj^Sj 

2 A Si 


,(* = 1 , 2 ) 


( 11 ) 


are the effective atom-photon coupling strengths. 

When choosing g$JA u = g 2 s JA Si , Q ri g r JA ri = 
fl Si g Si /A Si , and tp ai = -ip n = <pi, the Hamiltonian © 
turns into 


H — ujaQ^& ± c± ujqJ z ± Ai (a f + a) (12) 

x (j+e-^ 1 + J_e ivi ) + A 2 (J+e _iV2 + J_e^ 2 )(6 t + b ), 

where the effective atom-photon coupling strengths be¬ 
come 


1 9 si 

1 

2 A Sl 

” 2 A P1 

1 ^2 S 2 9s2 

_ 1 ^ 2 r 2 9r 2 

2 Ag 2 

~ 2 A ra 


(13) 


(14) 


If setting (fi = 0 and ip 2 = — 7 t/2 , the Hamiltonian (H 2 j) 
becomes 


H = uja^u ± oosb^b ± ujq J z (15) 

± Ai Ja,(a^ + cl) + A2 Jy(b^ + b). 

The Hamiltonian m is our required two-mode Dicke 
model, based on recent experiments of cavity-assisted Ra¬ 
man transitions PEI. In contrast to the convectional 
two-mode Dicke model achieved in the two-level atoms, 
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the Hamiltonian (fl5l) has a distinct property that all pa¬ 
rameters can be tuned independently. For example, the 
effective cavity frequencies uja and uo b depend on the de¬ 
tunings A a and A&, respectively; see Eqs. m and <nu- 
Thus, they can range from the positive to the negative. 
The choice of the different cavity frequencies uja and ujb 
help us to create a tunable two-axis spin Hamiltonian, as 
will be shown below. The effective atomic resonant fre¬ 
quency cjo can also be controlled by the detuning Ai; see 
Eq. ([8]). In addition, the effective atom-photon coupling 
strengths Ai and A 2 can be driven by the Rabi frequencies 
of Raman lasers; see Eqs. m and (HTj) . 


D. Generalized two-axis spin Hamiltonian 


In the following, we mainly consider the dispersive 
regime, i.e., {\wa\ , \^b\} {Ai, A 2 }. In such case, the 
photons are virtually excited, and we can use the Heisen¬ 
berg equations of motion [44j|, a = — i(Ai J x + ujao) = 0 
and b = —i(\ 2 J y + wsb) = 0, to obtain a = —Ai J x /ooa 
and b = — A 2 J v /ujb- As a result, the Hamiltonian (fl5l) 
becomes 


where 


Q(Jx + xjy) + W 0 Jz, 

(16) 

A? 


q =-, 

(17) 


u A \ 2 

(18) 

X ~ C 0 B Xj' 


In the Hamiltonian (fl6|) , the parameter q determines the 
nonlinear spin-spin interaction J 2 induced by the virtual 
photon, and the dimensionless parameter x reflects the 
ratio between the different nonlinear spin-spin interac¬ 
tions J 2 and J 2 . Due to existence of these nonlinear 
spin-spin interactions with the dimensionless parameter 
X and the effective atomic resonant frequency here 
we call the Hamiltonian ceu as a generalized two-axis 
spin Hamiltonian. In addition, equations ©, (Hzd, and 
(fl8|) show clearly that all the parameters, including the 
interaction strength g, the dimensionless parameter x> 
and the effective atomic resonant frequency ujq, can also 
be tuned independently in experiments. 

When chosen reasonable parameters, the generalized 
two-axis spin Hamiltonian m can reduce to some well- 
studied Hamiltonians. For example, when x > 0, the 
Hamiltonian ceu is a standard Lipkin-Meshkov-Glick 
model [45l-[47|| . When ujq = 0, the Hamiltonian (HU 
reduces to a generalized two-axis twisting Hamiltonian 
Hgtat = q{J 2 +xJy)- If further setting x = — 1 , a stan¬ 
dard two-axis twisting Hamiltonian Ht at = q(J 2 ~ Jy) 
is derived. Finally, when x = 0, the Hamiltonian (fl6|) 
turns into a generalized one-axis twisting Hamiltonian 
Hgoat = qJx + cjo J z 0, which reduces to the stan¬ 
dard one-axis twisting Hamiltonian -Ho at = qJ 2 f° r 



FIG. 3: (Color online). Numerical plots of the maximal 
squeezing factors as a function of the atomic number N 
for the different dimensionless parameters x* Insert: Numer¬ 
ical plot of the maximal squeezing factor as a function of 
the dimensionless parameter x> when the atomic number is 
chosen as N — 100. In both subfigures, the initial states are 
chosen as | J z = —j ). 


co>o = 0- These results implies that the Hamiltonian (fl6l) 
has an important application in achieving the required 
spin squeezing. 

Notice that when 00 a < 0, the effective spin-spin in¬ 
teraction strength q > 0. Thus, we can use the initial 
state | J z == — j) to discuss spin squeezing of the gener¬ 
alized two-axis spin Hamiltonian (H6j) . This initial state 
| J z = — j) can be easily prepared in experiments. In ad¬ 
dition, by combined with the effective atomic resonant 
frequency cjo 5 the dimensionless parameter x plays an 
important role in spin squeezing, as will be shown. 


III. SPIN SQUEEZING FACTOR 

In order to investigate spin squeezing, it is very neces¬ 
sary to consider the following time-dependent squeezing 
factor [H: 


£s (t) = min [AJlA*)] > (!9) 

where n± refers to an axis, which is perpendicular 
to the mean-spin direction no = J/ \J\ with \J\ = 

\J {Jx ) 2 + {Jy ) 2 + {Jz) 2 , and A A 2 = (A 2 ) - (A) 2 is the 

standard deviation. If |£|(t)| < 1 , the spin state is 
squeezed, and vice versa. 

In the spherical coordinates, no = 
(sin 6 cos (/?, sin 6 sin g?, cos 0 ) , where 0 = 

arccos((J^) / | J|), and p = arccos ((J x ) / | J\ sin#) for 
(J y ) > 0 or p = 2 tt — arccos (( J x ) / \J\ sin#) for (J y ) < 0. 
Two orthogonal bases are given by n\ — (— sin g?, cos g?, 0) 
and fi 2 = (— cos # cos g?, cos# sing?, — sin#). Thus, 

Jfii 7 2 = J ’ ^ 1,25 Jfi± = Jfti cos cf) T ? Jfi 2 sin0, and 
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FIG. 4: (Color online). Numerical plot of the maximal squeez¬ 
ing factor as a function of the effective atomic resonant 
frequency ujo. Insert: Numerical plot of the time-dependent 
squeezing factor £§(£). In both subfigures, the initial states, 
the dimensionless parameter, and the atomic number are cho¬ 
sen as | J z — — j), x — — 1, an d N = 100, respectively. 

min(Aj2^) could be achieved when p varies from 0 to 
27r in the plane that is perpendicular to the mean-spin 
direction no- It should be noticed that in experiments, 
the maximal squeezing factor ^ is usually measured 
0 0. Thus, in the following discussions, we mainly 
focus on this physical quantity. 

Before proceeding, we check the validity of the Hamil¬ 
tonian (USD, when the initial state is chosen as | J z = —j). 
For the generalized two-axis spin Hamiltonian, it is very 
hard to obtain analytical result of the spin squeezing fac¬ 
tor £§(£) 0. In Fig. [2j we numerically plot the corre¬ 
sponding spin squeezing factors (t) of the Hamiltio- 
nians m and m- It can be seen clearly that the 
results of the Hamiltonian (fT6|) are almost identical to 
those of the Hamiltonian (USD. Therefore, we will ap¬ 
ply the Hamiltonian (USD to discuss the experimentally- 
measurable maximal squeezing factor ^ in the rest of 
this paper. 


IV. MAXIMAL SQUEEZING FACTOR 

We first address a simple case without the effective 
atomic resonant frequency (cjo = 0): in which the gener¬ 
alized two-axis spin Hamiltonian (H6) ) reduces to the gen¬ 
eralized two-axis twisting Hamiltonian FFgtat = q(Jx + 
xJy)- I n Fig- 0 we numerically plot the maximal squeez¬ 
ing factor ^ of the Hamiltonian Hq tat as a function of 
the atomic number N for the different dimensionless pa¬ 
rameters x, when the initial state is chosen as | J z = — j). 
This figure shows clearly that when x = 0, the general¬ 
ized two-axis twisting Hamiltonian Hq tat becomes the 
standard one-axis twisting Hamiltonian iFoAT = qJxj 
whose maximal squeezing factor ^ scales as N~ 2 ' 3 
0. When increasing the dimensionless parameter x? the 



FIG. 5: (Color online). Numerical plots of the maximal 
squeezing factors as a function of the effective atomic reso¬ 
nant frequency ujq for the different dimensionless parameters 
X- The initial state and the atomic number are chosen as 
| J z = — j ) and N — 100, respectively. 

maximal squeezing factor decreases, i.e., spin squeez¬ 
ing is enhanced. In particular, when x = — 1, the Hamil¬ 
tonian Hq tat turns into the standard two-axis twisting 
Hamiltonian FFr at = q(~ Jy), whose maximal squeez¬ 
ing factor scales as V -1 [l|, as expected. In addition, 
the maximal squeezing factor ^ as a function of the di¬ 
mensionless parameter x is also plotted in the insert part 
of Fig. [3j This figure shows that x = — 1 is an optimal 
point to achieve the maximal squeezing factor of the 
generalized two-axis twisting Hamiltonian FFgtat- 

In real experiments, the effective atomic resonant fre¬ 
quency cjo always exists. In Fig. 0J we numerically plot 
the maximal squeezing factor of the Hamiltonian 
H = q(J 2 — J 2 ) + ujqJ z as a function of the effective 
atomic resonant frequency cjo, when the the initial state 
is chosen as \J Z = — j). It can be seen from this figure 
that with the increasing of the effective atomic resonant 
frequency cjo in the standard two-axis twisting Hamilto¬ 
nian FFr at? the maximal squeezing factor ^ increases, 
i.e., spin squeezing is reduced. 

From the above discussions, we argue that when in¬ 
creasing the dimensionless parameter x (from x = — 1) 
or introducing the effective atomic resonant frequency ujq 
in the standard two-axis twisting Hamiltonian FFr at, the 
maximal squeezing factor increases, i.e., spin squeez¬ 
ing is reduced. Surprisingly, when we control these two 
parameters simultaneously, spin squeezing can be en¬ 
hanced largely. To see this clearly, in Fig. [5] we nu¬ 
merically plot the maximal squeezing factors ^ of the 
generalized two-axis spin Hamiltonian ([T 6]h i.e., H = 
q(J 2 + xJy) + <^o Jzi as a function of the effective atomic 
resonant frequency ujq for the different dimensionless pa¬ 
rameters x, when the initial state is chosen as | J z = — j). 
This figure shows that in the case of x = —0.05 or 
X = —0.5, when increasing the effective atomic resonant 
frequency cjq, the maximal squeezing factor ^ first de- 
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creases largely and then increases, i.e., spin squeezing is 
first enhanced largely and then reduced. For the gen¬ 
eralized one-axis twisting model Hq oat = + woJ z , 

the maximal squeezing factor ^ has a similar behavior 
(see green dash-dot line in Fig. © 0, but its magnitude 
cannot arrive at the order of our considered two-axis spin 
Hamiltonian (fl6l) because these two Hamiltonians have 
different scalings with respect to the atomic number N. 
These results are benefit for achieving the required spin 
squeezing in experiments. 

V. DISCUSSIONS 

Finally, we estimate the relative parameters, based on 
recent experiments. Two stable ground states are cho¬ 
sen respectively as |Go) = \F = l,rnp = 1 ) and |Gi) = 
| F = 2 ,mi? = 2 ) of ultracold rubidium 87 atoms, where 
F is the total angular momentum and is the mag¬ 
netic quantum number. It means that the atomic decay 
rate is 7 / 2tt = 3.0 MHz, which is the same order as the 
photonic decay rate k,/2tt = 1.3 MHz [4l], [42]. In such 
system, the atom-photon coupling strengths can reach 
g Sl /27r = 20 MHz and g S2 /2 tt = 15 MHz, respectively. 
When CI Si /2tt = 20 MHz and A Sl /2i r = 100 MHz, which 
is responsible for deriving Eq. (flTj) . \\j2i r = 2 MHz, 
and thus q/2 tt = 0.2 MHz for da/ 27t = —20 MHz. This 
choice of the detuning uja also satisfies the dispersive con¬ 
dition, which is a key condition to realize our generalized 
two-axis spin Hamiltonian (H 6 j) . For the dimensionless 
parameter y, it can be easily controlled by tuning both 
the detuning ujb and the Rabi frequency U S2 . When 
N = 100, numerical result shows that the shortest time 
for generating the maximal squeezing factor ^ is about 
t m ~ 20 ns [see, for example, in the insert part of Fig.0], 
which is shorter than both the atomic and photonic life¬ 


times 7 -1 and ft -1 . With the increasing of the atomic 
number N and the atom-photon coupling strength g Sl , 
£ m becomes shorter and shorter. This indicates that our 
proposal can be accessible in the current experimental 
setups. 


VI. CONCLUSIONS 

In summary, we have proposed an experimentally- 
feasible system, in which an ensemble of ultracold six- 
level atoms interacts with two quantized cavity fields and 
two pairs of Raman lasers, to realize a generalized two- 
axis spin Hamiltonian H = q(J% + xJy) + woJ z . We have 
numerically calculated the experimentally-measurable 
maximal squeezing factor and revealed that when ujq = 0 
and % = — 1 , the maximal squeezing factor ^ scales as 
V -1 . More importantly, we have found that by com¬ 
bined with the dimensionless parameter x(> — 1 ), the 
effective atomic resonant frequency uoq can enhance spin 
squeezing largely. Our results are benefit for achieving 
the required spin squeezing in experiments, and have a 
potential application in quantum information and quan¬ 
tum metrology. 
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